Neutron- 3 H scattering above the four-nucleon breakup threshold 
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The four-body equations of Alt, Grassberger and Sandhas are solved for the neutron- 3 H scattering 
at energies above the four-nucleon breakup threshold. The accuracy and practical applicability of 
the employed complex energy method is significantly improved by the use of integration with the 
special weights. This allows to obtain fully converged results with realistic nuclear interactions. A 
satisfactory description of the existing neutron- 3 H elastic scattering data is obtained. 
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The four-nucleon reactions is an ideal but also highly 
challenging field to test few-nucleon interaction models. 
The problem of elastic nucleon-trinucleon scattering be- 
low the inelastic threshold has already been solved with 
high accuracy using several ab initio methods with realis- 
tic nuclear potentials. These methods include the hyper- 
spherical harmonics (HH) expansion [H-Q; the Faddeev- 
Yakubovsky (FY) equations Q for the wave function 
components in the coordinate space 0, @|, and the Alt, 
Grassberger and Sandhas (AGS) equations H, M for the 
transition operators in the momentum space [9h11| . A re- 
cent benchmark 12j reported a good agreement between 
the HH, FY, and AGS techniques for the neutron- 3 H 
(n- 3 H) and proton- 3 He (p- 3 He) scattering. Furthermore, 
deuteron-deuteron (d-d) collisions, including the transfer 
reactions to p- 3 H and n- 3 He final states, have been cal- 
culated using the resonatin g-gr oup method (RGM) [l3[ 
and the AGS framework [lj,ll5|. However, also these cal- 
culations were limited to energies below the three-cluster 
breakup threshold. At higher energies, especially above 
the four-body breakup threshold, the asymptotic bound- 
ary conditions in the coordinate space become nontriv- 
ial due to open two-, three- and four-cluster channels. 
In the momentum-space framework one is faced with a 
very complicated structure of singularities in the kernel 
of integral equations. Formally, these difficulties can be 
avoided by rotation to complex coordinates [16| or con- 
tinuation to complex energy [l?], EH that lead to bound- 
state like boundary conditions and nonsingular kernels. 
However, technical complications may arise in practical 
calculations. Indeed, the applications to the four-nucleon 
scattering so far have been very limited [19, 20] and none 
of them uses realistic interactions. On the other hand, 
the no-core shell model RGM [2l[ includes in the model 
space only the ground state of the three-nucleon system 
which is insufficient. In Ref. [22[ this shortcoming was 
partially corrected by adjusting the predictions to the 
experimental data. 

The aim of the present work is to overcome the 
above limitations by performing realistic well-converged 
four-nucleon scattering calculations above the four-body 
breakup threshold. We use the complex energy method 
[l7| but introduce important technical improvements. 
Although in the AGS framework employed by us the 
Coulomb force can be included via the screening and 



renormalization method [23|, [24| , the present numerical 
results are restricted to the Coulomb-free n- 3 H case. 

We treat the nuclcons as identical particles in the 
isospin formalism and therefore use the AGS equations 
for the symmetrized four-particle transition operators 
Uf3 a as derived in Ref. (§), i.e., 

U\\ = — (Got Go) 1 ^34 — P34U1G0 1 Go U11 

+ U 2 G tG U 21 , (la) 

W 2 i = (G tGo)- 1 (l-P 3 4) + (l-P34)C/iGotG W 11 , 

(lb) 

^12 = (Go t Go) 1 — P34C/1G0 1 Go U12 + U2G0 1 GqJA-22, 

(lc) 

U 2 2 = {l-P3i)U l G tG U l2 . (Id) 

Here, a — 1 corresponds to the 3 + 1 partition (12,3)4 
whereas a = 2 corresponds to the 2+2 partition (12)(34); 
there are no other distinct two-cluster partitions in the 
system of four identical particles. 



G = {Z- Ho)- 1 



(2) 



is the free resolvent with the complex energy parameter 
Z = E + is and the free Hamiltonian Hq, 



t = V + vGnt 



(3) 



is the pair (12) transition matrix derived from the poten- 
tial v, and 



U a — P a G 1 + P a t Gq U a , 



(4) 



are the symmetrized 3+1 or 2+2 subsystem transition 
operators. The basis states are antisymmetric under ex- 
change of two particles in the subsystem (12) for the 3 + 1 
partition and in (12) and (34) for the 2 + 2 partition. The 
full antisymmetry of the four-nucleon system is ensured 
by the permutation operators P a b of particles a and b 
with P x = P 12 P23 + Pw P23 and P 2 = P13 P 24 . 

The scattering amplitudes for two-cluster reactions at 
available energy E = e a + p^,/2fi a — ep +p|/2/z^ are 
obtained as the on-shell matrix elements (pp\T/3 a \p a ) — 
<S/3a(<Aa|^9a|<^a) in the limit e — > +0. Here \(p a ) is the 
Faddeev component of the asymptotic two-cluster state 
in the channel a, characterized by the bound state energy 
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e a < 0, the relative momentum p Q , and the reduced mass 
/i Q . Thus, depending on the isospin, e\ is the ground 
state energy of 3 He or 3 H, and €2 is twice the deuteron 
energy e d - Sp a are the symmetrization factors @, e.g., 
5n = 3. The amplitudes for breakup reactions are given 
by the integrals involving Ufia\<t>a) [2a |. 

We solve the AGS equations (fTJ in the momentum- 
space partial-wave framework. The states of the total 
angular momentum J with the projection M. are 



k X ky k z ^l z ( K \ly\{l x S X )j X . 



defined as 
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2 + 2. Here k x , k y and k z are the four-particle Jacobi 
momenta in the convention of Ref. [2^], l x , l y , and l z are 
the associated orbital angular momenta, j x and j y are 
the total angular momenta of pairs (12) and (34), J y is 
the total angular momentum of the (123) subsystem, s y 
and s z are the spins of nucleons 3 and 4, and S x , S y , and 
S z are channel spins of two-, three-, and four-particle 
system. A similar coupling scheme is used for the 
isospin. In the following we abbreviate all discrete 
quantum numbers by v. The reduced masses associated 
with Jacobi momenta k x and k y in the partition a will 
be denoted by \i ax and fi ay , respectively. 

An explicit form of integral equations is obtained by 
inserting the respective completeness relations 
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(5) 

between all operators in Eqs. (JTJ. The integrals are 
discrctized using Gaussian quadrature rules 26] turn- 
ing Eqs. ((J) into a system of linear equations as de- 
scribed in Ref. jfj}. However, in the limit e — > +0 
needed for the calculation of the observables the kernel 
of the AGS equations contains integrable singularities. 
At E + ie — e a — k\j2[i a — > the subsystem transition 
operator in the bound state channel has the pole 
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Furthermore, at E 



Pa 1 0a ) S aa (0a | Pa 

E + ie — e a — k z /2fj, a 



(6) 



ie - e d - ky/2fj, ay - k z /2fj, a -> 
the two-nucleon transition matrix in the channel with 
the deuteron quantum numbers for the pair (12) has the 
pole 



v\(/>d)(<i>d\v 



E + ie - e d - kl/2^jL ay - k 2 z /2^ a ' 



(7) 



with \<j>d) being the pair (12) deuteron wave function. 
Finally, the free resolvent (|2J) obviously becomes singular 
at E + ie - k 2 x j2\i ax - k y /2fi ay - k 2 z /2^i a -t 0. 

At energies below the three-cluster threshold only © 
singularities are present that in our previous calculations 
Q were treated reliably by the subtraction technique. 
However, above the four-body breakup threshold all three 
kinds of singularities are present. Their interplay with 
permutation operators and basis transformations leads 



to a very complicated singularity structure of the AGS 
equations. As proposed in Ref. [17J, this difficulty can 
be avoided by performing calculations for a set of finite 
e > values where the kernel contains no singularities 
and then extrapolating the results to the e — > +0 limit. 
However, this extrapolation is only precise with not too 
large e values. On the other hand, for small e the kernel 
of the AGS equations, although formally being nonsingu- 
lar, may exhibit a quasi-singular behavior thereby requir- 
ing dense grids for the numerical integration. This is no 
problem in simple model calculations with rank-one sep- 
arable potential and very few channels pTj| where one can 
use a large number of grid points. However, in practical 
calculations with realistic potentials and large number 
of partial waves necessary for the convergence one has to 
keep the number of integration grid points possibly small 
and therefore a more sophisticated integration method is 
needed. 

We take over from Refs. [13) EH the idea of the com- 
plex energy method and the e —> +0 extrapolation pro- 
cedure (analytic continuation via continued fraction) but 
we introduce an important technical improvement when 
calculating Ufj a at finite e. We use the method of special 
weights for numerical integrations involving any of the 
above-mentioned quasi-singularities, i.e., 
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-dxfn^2f{xj)wj(n,x ,yo,a,b). (8) 
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The quasi-singular factor (xg + iyo — x 11 )" 1 is sepa- 
rated and absorbed into the special integration weights 
Wj(n, xo, yo, a, b). The set of TV grid points {xj} where 
the remaining smooth function f(x) has to be evalu- 
ated is chosen the same as for the standard Gaussian 
quadrature. However, while the standard weights are 
real [26j], the special ones Wj(n,xo,yo,a,b) are complex. 
They are chosen such that for a set of N test functions 
fj(x) the result (0 is exact. A convenient and reliable 
choice of {fj(x)} are the N spline functions {Sj(x)} re- 
ferring to the grid their construction and properties 
are described in Refs. j26H28|. The corresponding special 
weights are 



Wj(n,x ,yo,a,b) 



Sj(x) 



Xq + Wo ~ x n 



-dx, 



(9) 



where the integration can be performed either analyti- 
cally or numerically using sufficiently dense grid. This 
choice of special weights guarantees accurate results for 
quasi-singular integrals ([SJ with any f(x) that can be ac- 
curately approximated by the spline functions {Sj(x)}. 

In the integrals over the momentum variables one has 
n = 2, a = 0, and b — > 00. For example, when 
solving the Lippmann-Schwinger equation the inte- 
gration variable in Eq. @ is the momentum k x with 
xl = 2[i ax (E - k 2 /2^i ay - k\j2^ a ) and y = 2^ ax e. 
Alternatively, the quasi-singularity can be isolated in a 
narrower interval < a < b < 00 and treated by special 
weights only there. 
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n( 3 Po) 


[1.0,2.0] 


62.63 


0.990 


43.03 


0.959 


65.27 


0.950 


[1.2,2.0] 


62.60 


0.991 


43.04 


0.959 


65.29 


0.951 


[1.4,2.0] 


62.67 


0.991 


43.03 


0.958 


65.27 


0.950 


[1.2,1.8] 


62.65 


0.992 


43.03 


0.959 


65.28 


0.950 


1.4 


73.37 


0.916 


44.77 


0.840 


67.38 


0.933 



TABLE I. Elastic phase shifts (in degrees) and inelastici- 
ties in selected partial waves for n- 3 H scattering at 22.1 MeV 
neutron energy. Results for INOY04 potential obtained using 
different sets of e values ranging from e m i n to e max (in MeV) 
are compared. In the last line the predictions with e = 1.4 
MeV without extrapolation are given. 



Other numerical techniques for solving the four- 
nucleon AGS equations are taken over from Ref. 
They include Pade summation 29] of Neumann series 
for the transition operators U a and Up a using the al- 
gorithm of Ref. [30| and the treatment of permutation 
operators (basis transformations) using the spline inter- 
polation. The specific form of the permutation operators 
P leads to a second kind of quasi- singular integrals © 
with n = 1, a = —1, b = 1, and the integration variable 
x = Ity ■ k y or k' z ■ k z being the cosine of the angle between 
the respective initial and final momenta. 

We note that the above integration method is not suf- 
ficient in the vanishing e limit since for n = 1 and ya = 
the result of the integral © contains the contribution 
f(xo) ln[(i£o + l)/(^o — 1)] with logarithmic singularities 
at xq = ±1. At finite small e the result of ([5]) may ex- 
hibit a quasi-singular behavior. However, since the loga- 
rithmic quasi-singularity is considerably weaker than the 
pole quasi-singularity, at not too small e it is sufficient 
to use the standard integration. 

We start by applying the complex energy method with 
special integration weights to the n- 3 H scattering below 
the three-cluster threshold where our previous results [9| 
obtained at real energies using subtraction technique are 
available for comparison. In the test calculations at 3.5 
and 6.0 MeV neutron energy we find a very good agree- 
ment between the two methods, better than 0.05% for 
all relevant phase shifts and observables. The considered 
e values range between 0.2 and 2.0 MeV, and the num- 
ber of grid points is not increased as compared to the 
real-energy calculations @- 

Next we test the numerical reliability of our technique 
above the four-nucleon breakup threshold. We use a real- 
istic dynamics, namely, the high-precision inside-nonlocal 
outside- Yukawa (INOY04) two-nucleon potential by Do- 
leschall 0, |3l| that reproduces experimental binding en- 
ergies of 3 H (8.48 MeV) and 3 He (7.72 MeV) without 
an irreducible three-nucleon force. We consider a large 
number of four-nucleon partial waves sufficient for the 
convergence, i.e., l x ,ly,lz,jx,jy,J y < 4 and J < 5. In- 
cluding more partial waves yields only entirely insignif- 
icant changes. There are too many numerical param- 
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FIG. 1. (Color online) Differential cross section and neu- 
tron analyzing power for elastic n- 3 H scattering at 22.1 MeV 
neutron energy as functions of cm. scattering angle. Re- 
sults obtained using different sets of e values ranging from 
e m in to £ max with the step of 0.2 MeV are compared; they 
are indistinguishable. The dotted curves refer to the e = 1.4 
MeV calculations without extrapolation that have no physical 
meaning but show the importance of the extrapolation. 

eters (numbers of points for various integration grids) 
to demonstrate the stability of our calculations with re- 
spect to each of them separately. We found that 10 grid 
points are sufficient for all angular integrations but 30 
to 40 grid points are needed for the discretization of Ja- 
cobi momenta. The e — > +0 extrapolation yields stable 
results only if sufficiently small e are considered and at 
each of them the respective calculations are numerically 
well converged. We therefore study in Fig. Q] the sta- 
bility of the e —> +0 results obtained via extrapolation 
using different e sets ranging from e m j n to e max with the 
step of 0.2 MeV. We show the differential cross section 
da/dQ, and neutron analyzing power A y for elastic n- 3 H 
scattering at E n = 22.1 MeV neutron energy. We find a 
very good agreement between the results obtained with 
[1.0,2.0], [1.2,2.0], [1.4,2.0], and [1.2,1.8] 
MeV, confirming the reliability of our calculations. In 
addition, we show in Fig. [T] the predictions referring to 
e = 1.4 MeV without extrapolation that don't have phys- 
ical meaning. The difference between e — > +0 and e = 1.4 
MeV results demonstrates the importance of the extrapo- 
lation. Furthermore, in TableUwe collect the correspond- 
ing values for selected phase shifts <5 and inelasticities r], 
i.e., we parametrize the elastic S'-matrix as s = rje 2lS . 
As already can be expected from Fig. [TJ the stability of 
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FIG. 2. (Color online) Differential cross section for elastic 
n- 3 H scattering at 14.1, 18.0, and 22.1 MeV neutron energy. 
Results obtained with INOY04 (solid curves) and CD Bonn 
(dashed-dotted curves) potentials are compared with the ex- 
perimental data from Refs. 22, 32, 33]. 



the results with respect to changes in [e m i n , e ma x] is very 
good. The variations are slightly larger in the 1 5o state 
where also the difference between the finite e and e — > +0 
results is most sizable. Nevertheless, from Fig. [1] and Ta- 
ble U one can conclude that with a proper e choice as 
few as four different e values are sufficient to obtain the 
physical e — > +0 results with good accuracy. 

For curiosity, in J = states we performed the cal- 
culations keeping the same grids but with standard in- 
tegration weights. We found that they fail completely 
at e values from Table HI with the errors of the e — > +0 
extrapolation being up to 10 % for phase shifts and up 
to 25 % for inelasticity parameters. On the other hand, 
at large e > 4 MeV the two integration methods agree 
well but the e —> +0 extrapolation has at least one order 
of magnitude larger inaccuracies than those in Table HI 

After establishing the reliability of our calculations we 
proceed to the comparison with the experimental data. 
In addition to the INOY04 potential we present results 
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FIG. 3. (Color online) Neutron analyzing power for elastic 
n- H scattering. INOY04 and CD Bonn predictions at E n — 
22.1 MeV are compared with the data from Ref. [H. INOY04 
results at 14.1 MeV and 18.0 MeV are also shown. 



derived from the CD Bonn potential [3J] that underbinds 
the 3 H nucleus by 0.48 MeV. In Fig. [H we show the dif- 
ferential cross section for elastic neutron- 3 H scattering at 
14.1, 18.0, and 22.1 MeV neutron energy. Except for the 
minimum around 115 degrees, the predictions are insen- 
sitive to the choice of the potential. At E n = 14.1 MeV 
the new data set by Frenje et al. [22j is described very 
well. Other existing data at this energy are consistent 
with Ref. (22| but have larger error bars; we only show 
the data by Debertin et al. [H|. At 18.0 MeV the data 
sets by Debertin et al. [32[ and Seagrave et al. [33| are 
inconsistent with each other around the minimum while 
the theoretical predictions lie in the middle. The results 
at E n — 22.1 MeV are compared with the data taken 
at 21 and 23 MeV by Seagrave et al. (33[. The predic- 
tions lie between the two data sets except for the mini- 
mum region. However, given the agreement between the 
22 1 and [HJ data and disagreement between the [HI and 
32| data, one may question the reliability of the data 
by Seagrave et al. in the minimum region. Thus, new 
measurements are needed to resolve this discrepancy. 
In Fig. [3] we present the neutron analyzing power for 



elastic n- H scattering at E n — 22.1 MeV. To study the 
energy dependence we also show INOY04 predictions at 
E n = 14.1 and 18.0 MeV. We observe that the sensitiv- 
ity to the nuclear force model and energy is considerably 
weaker as c omp ared to the regime below three-cluster 
threshold Most remarkably, in contrast to low en- 

ergies where the famous p- 3 He A y -puzzle exists fll llol . l35j . 
the peak of A y around 120 degrees is described very well 
but there is a discrepancy in the minimum region. This 
is somehow similar to the three-nucleon system where 
the nucleon-deuteron A y -puzz\e existing at low energies 
disappears as the energy increases (36j . 

In this work we do not calculate explicitly the breakup 
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at Ref. 


14.1 


928 


19 


947 


913 


28 


941 


978 ±70 [3JJ 


18.0 


697 


41 


738 


689 


48 


737 


750 ±40 [37] 


22.1 


536 


61 


597 


524 


70 


594 


620 ±24 [38] 



TABLE II. n- 3 H elastic <r e , breakup ab, and total <7t cross 
sections (in mb) at selected neutron energies (in MeV). 

reactions. However, the total n- 3 H cross section <r t = 
<r e + <j b with the elastic <r e and three- and four-cluster 
breakup a b contributions is calculated using the optical 
theorem. The results at three considered energies are 
collected in Table [TT1 The at predictions are slightly lower 
than the experimental data from Refs. (37l . [38[ but in 
most cases they agree within error bars. The breakup 
cross section a b is most sensitive to the potential model, 
probably due to differences in 3 H binding energy and 



breakup threshold positions. 

In summary, we performed fully converged neutron- 3 H 
scattering calculations with realistic potentials above the 
four-nucleon breakup threshold. The symmetrized Alt, 
Grassberger, and Sandhas four-particle equations were 
solved in the momentum-space framework. We greatly 
improved the accuracy and the efficiency of the complex 
energy method by using numerical integration technique 
with special weights for all quasi-singularities of the pole 
type. The calculated elastic scattering observables show 
little sensitivity to the interaction model. The overall 
agreement with the experimental data for the elastic dif- 
ferential cross section and neutron analyzing power as 
well as for the total cross section is quite good, consider- 
ably better than in the low-energy n- 3 H and p- 3 He scat- 
tering. Extension of the method to other reactions in the 
four-nucleon system is in progress. 

The author thanks R. Lazauskas and A. C. Fonseca 
for valuable discussions and J. A. Frenje for providing 
the data. 
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